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MACRO constraints on violation of Lorentz invariance 
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The energy spectrum of neutrino-induced upward-going muons in MACRO has been analysed in terms of 
relativity principles violating effects, keeping standard mass-induced atmospheric neutrino oscillations as the 
dominant source of — > v T transitions. The data disfavor these exotic possibilities even at a sub-dominant level, 
and stringent 90% C.L. limits are placed on the Lorentz invariance violation parameter | Ad| < 6 x 10 -24 at sin 26 v 
= and |Au| < 2.5 -f- 5 x 1(T 26 at sin 26^ = ±1. These limits can also be re-interpreted as upper bounds on the 
parameters describing violation of the Equivalence Principle. 



1. Introduction 

Neutrino mass-induced oscillations are the best 
explanation of the atmospheric neutrino problem 
|l|2|3|4j . Two flavor — > v T oscillations are 
strongly favored over a wide range of alterna- 
tive solutions [5|6|l7[8|9|10|llj . These alternative 
mechanisms have been considered under the hy- 
pothesis that each one of them solely accounts for 
the observed effects. We address the possibility 
of a mixed scenario: one mechanism, the mass- 
induced flavor oscillations, is considered domi- 
nant and a second mechanism is included in com- 
petition with the former. We studied, as sub- 
dominant mechanism, neutrino flavor transitions 
induced by violations of relativity principles, i.e. 
violation of the Lorentz invariance (VLI) or of the 
equivalence principle (VEP). 

In this mixed scenario, we assume that neu- 
trinos can be described in terms of three dis- 
tinct bases: flavor eigenstates, mass eigenstates 
and velocity eigenstates, the latter being charac- 
terized by different maximum attainable veloc- 
ities (MAVs), and consider that only two fam- 
ilies contribute to the atmospheric neutrino os- 
cillations. When both mass-induced and VLI- 
induced neutrino oscillations are considered si- 
multaneously, the survival probability can be 
expressed as [9110] 

sin 2 26 sin 2 



where and Q are given by 
20 = arctan(ai/<22); il 
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The terms a\ and ct2 contain the relevant phys- 
ical information and the whole domain of vari- 
ability of the parameters can be accessed with 
the requirements Am 2 > 0, < 9 m < tt/2, 
Av > and -tt/4 < 6 V < it /4. The func- 
tional form of the oscillation probabilities for 
mass induced (VLI induced) oscillations exhibits 
an L/ E v (LE U ) dependence on the neutrino ener- 
gies and pathlengths respectively. The same for- 
malism also applies to violation of the equivalence 
principle, after substituting Aw/2 with the adi- 
mensional product |</>|A7; A7 is the difference of 
the coupling constants for neutrinos of different 
types to the gravitational potential <f> |12j . 

As shown in [TO], the most sensitive tests of 
VLI can be made by analysing the high energy 
tail of atmospheric neutrinos at large pathlength 
values. As an example, Fig. [1] shows the energy 



dependence of the v„ 



survival probability 



as a function of the neutrino energy, for neutrino 
mass-induced oscillations alone and for both mass 
and VLI-induced oscillations for Aw = 2 x 10 -25 
and sin 2 29 v = ±1. Notice that for m v < 1 
eV, neutrinos with energies larger than 100 GeV 
are extremely relativistic, with Lorentz 7 factors 
larger than 10 11 . 



2. Experimental data and analysis 

In order to analyse the MACRO data in terms 
of VLI, we used a subsample of 300 upward- 
throughgoing muons whose energies were esti- 
mated via multiple Coulomb scattering in the 7 
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|Av = 2 10 2b , sin2 9 v =±T 




Figure 1. Energy dependence of the — ► 
survival probability for mass-induced oscillations 
alone (continuous line) and mass-induced + VLI 
oscillations for Av — 2 • 10~ 25 and sin 28 v = 
±1, dashed lines for positive values, dotted lines 
for negative values. The neutrino pathlength was 
fixed at L = 10000 km. 



horizontal rock absorbers in the lower appara- 
tus [13114] . The energy estimate was obtained 
using the streamer tubes in drift mode, which 
allowed to considerably improve the spatial res- 
olution of the detector (~3 mm). The overall 
neutrino energy resolution was of the order of 
100%, mainly dominated by muon energy losses 
in the rock below the detector (note that (E^) ~ 
0.4 (E v )). Upgoing muon neutrinos from this 
sample are particularly suited for sensitive tests 
of VLI and VEP, since they have large energies 
((E) ~ 50 GeV) and large zenith angles (> 120°), 
which correspond to a median value of neutrino 
pathlengths of about 10000 km. Following the 
analysis in Ref . [14] , we selected a low and a high 
energy sample by requiring that the reconstructed 
neutrino energy E r v ec should be E r v ec < 30 GeV 
and E r v ec > 130 GeV. The number of events sur- 
viving these cuts is Ni ow = 49 and Nhigh — 58, 
respectively; their median energies, estimated via 
Monte Carlo, are 13 GeV and 204 GeV (assum- 
ing mass- induced oscillations). The analysis then 
proceeds by fixing the neutrino mass oscillation 



parameters at the values of [2] . Then, we scanned 
the plane of the two free parameters (Av, V ) min- 
imizing the x 2 function comprehensive of statis- 
tical and systematic uncertainties [TS]. For the 
Monte Carlo simulation described in [14] the neu- 
trino fluxes in input is given by [16] . The largest 
relative difference of the extreme values of the 
MC expected ratio Ni ow /Nhi g h is 13%. How- 
ever, in the evaluation of the systematic error, 
the main sources of uncertainties for this ratio 
(namely the primary cosmic ray spectral index 
and neutrino cross sections) have been separately 
estimated and their effects added in quadrature 
(see [14] for details): in this work, we use a con- 
servative 16% theoretical systematic error on the 
ratio Ni ow /Nhi g h- The experimental systematic 
error on the ratio was estimated to be 6%. The 
inclusion of the VLI effect does not improve the 
X 2 in any point of the (Av, V ) plane, compared to 
mass-induced oscillations stand-alone, and proper 
upper limits on VLI parameters were obtained. 
The 90% C.L. limits on Aw and V) computed 
with the Feldman and Cousins prescription [17j . 
are shown by the dashed line in Fig. [21 The energy 
cuts described above (the same used in Ref. |14|L 
were optimized for mass-induced neutrino oscil- 
lations. In order to maximize the sensitivity of 
the analysis for VLI induced oscillations, we per- 
formed a blind analysis, based only on Monte 
Carlo events, to determine the energy cuts which 
yield the best performances. The results of this 
study suggest the cuts E r v ec < 28 GeV and E r v ec > 
142 GeV; with these cuts the number of events in 
the real data are N' low = 44 events and N' high = 
35 events. The limits obtained with this selection 
are shown in Fig. [3 by the continuous line. As 
expected, the limits are now more stringent than 
for the previous choice. 

An independent and complementary analysis 
was performed on a sample of events with a re- 
constructed neutrino energy 25 GeV < E r v ec < 75 
GeV. The number of events satisfying this condi- 
tion is 106. A negative log-likelihood function was 
built event by event and then fitted to the data. 
We allowed mass-induced oscillation parameters 
to vary inside the MACRO 90% C.L. region and 
we left VLI parameters free in the whole (Av, 9 V ) 
plane. The upper limit on the Av parameter re- 
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suiting from this analysis is slowly varying with 
Am 2 and is of the order of « 1CP 25 . 

3. Conclusions 

We have searched for "exotic" contributions 
to standard mass-induced atmospheric neutrino 
oscillations arising from a possible violation of 
Lorcntz invariance. Two different and comple- 
mentary analyses were performed on the data, 
both of them yielding compatible upper limits for 
the VLI contribution. We used a subsample of 
MACRO muon events for which an energy mea- 
surement was made via multiple Coulomb scatter- 
ing. The inclusion of VLI effects does not improve 
the fit to the data, and we conclude that these 
effects are disfavored even at the sub-dominant 
level [TS]- The VLI parameter bound is (at 90% 
C.L.) | At; | < 3 x 1CT 25 . This result may be rein- 
terpreted in terms of 90% C.L. limits of param- 
eters connected with violation of the equivalence 
principle, giving the limit |0A7| < 1.5 x 10~ 25 . 
The second approach exploits the information 
contained in a data sub-set characterized by inter- 
mediate muon energies. It is based on the maxi- 
mum likelihood technique, and considers the mass 
neutrino oscillation parameters inside the 90% 
border of the global result [5] . The obtained 90% 
CL upper limit on the Av VLI parameter versus 
the assumed Am 2 values is also around 10~ 25 . 
The limits reported in this paper are comparable 
to those estimated using Super-Kamiokande and 
K2K data [9]. 
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